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ABSTRACT 

We have carried out an ALMA survey of 15 confirmed or candidate low-mass (< 0.2 M Q ) members of the TW Hya Association (TWA) 
with the goal of detecting molecular gas in the form of CO emission, as well as providing constraints on continuum emission due to 
cold dust. Our targets have spectral types of M4-L0 and hence represent the extreme low end of the TWA’s mass function. Our ALMA 
survey has yielded detections of 1.3mm continuum emission around 4 systems (TWA 30B, 32, 33, & 34), suggesting the presence 
of cold dust grains. All continuum sources are unresolved. TWA 34 further shows 12 CO(2-l) emission whose velocity structure is 
indicative of Keplerian rotation. Among the sample of known —7—10 Myr-old star/disk systems, TWA 34, which lies just ~50 pc from 
Earth, is the lowest mass star thus far identified as harboring cold molecular gas in an orbiting disk. 
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1. Introduction 

Beginning with the identification of the TW Hya Association 
(TWA; Kastner et a l,|1997 1 several stellar associations with ages 
-8-200 Myr have been identified in close proximity to the Earth 
( |Zuckerman & Song||2004| |Torres et al.||2008| ). These young 
moving groups serve as excellent laboratories to explore the evo¬ 
lution of stellar and planetary properties (e.g., Marois et al. 2008[ 
Lagrange et al.|2010[ [Rodriguez et al.|2010] l. Of particular inter¬ 
est is the growing number of M-dwarfs identified as members 
of these groups and how their properties, such as disk lifetimes, 
compare with those of higher mass members. 

Surveys of star forming regions have shown that protoplan¬ 
etary disk lifetimes are on the order of -2-3 Myr (Williams & 
Cieza j20Tl| . To date, we know only of a handful of pre-main 
sequence stars within 100 pc of Earth that, despite having ages 
of -5-20 Myr, still host disks that appear primordial in nature. 
These include TW Hya (-8 Myr, -50 pc; |Kastner et al.|1997), 
V4046 Sgr (-20 Myr, -70 pc; |Kastner et al.||2008), MP Mus 
(-5-7 My r, -100 pc;|Kastner et al.|2010| >, and T Cha (-5-7 Myr, 
-110 pc; |Sacco et al.|2014[ >. All of these are roughly solar-mass 
(K type) stars; whether or not lower-mass (M type) stars can re¬ 
tain primordial disks to such relatively advanced ages remains to 
be determined. 

Young M-dwarfs can exhibit strong levels of UV and X-ray 
radiation, either as a result of chromospheric and coronal activity 
or active accretion. This high energy radiation can drive disk dis¬ 


sipation via photoevaporation (e.g., Gorti & Hollenbach 2009| . 
Nevertheless, the limited studies that have been performed to ex¬ 
plore the disk lifetimes of substellar objects have found disk dis¬ 
sipation timescales are at least as long as for solar-mass stars 
(see, e.g., |Luhman & Mamajek||2012| |Ercolano et al.||2011| 
Williams & Cieza|2011 and references therein). Investigations 


aimed at establishing the masses of gas and dust around similarly 
young M stars are key to advance our understanding of disk evo¬ 
lution timescales and processes. 

Among the nearby young moving groups, the TWA is par¬ 
ticularly interesting, as it represents an important evolutionary 
stage in the lives of protoplanetary disks. With an age of -8 Myr 
(Torres et al.|[2008[ [Ducourant et al .(2014| , the TWA represents 
an epoch coinciding with the time required for giant planet for¬ 
mation via core accretion (eg, |Chabrier et al.|2014[ >. The disks in 
the TWA range from evidently gas-poor debris disks to at least 
one example of a long-lived, gas-rich, apparently primordial disk 
— i.e., the disk orbiting TW Hya i tself (e.g.,|Riviere-M arichalar 
|et al.|2013]|Schneider et al.|2012aH Andrews et al.|2010t and ref- 

erences therein). Recent work has yielded many new M-dwarf 
members and candidate members of the T WA (|Looper et aL 
2007 2010a|b[ |Looper||201 1 [ |Shkolnik et ~aL)|2011 [ Rodriguez 


et al. 2011 1 Schneider et al.||2012b| Gagne et al. |2014b| >. With 

masses lower than -0.2 M a , these new mid/late M-dwarfs con¬ 
stitute a sample of objects that allows us to probe disks among 
the poorly explored substellar population, and to do so at the 
relatively advanced age of the TWA. 
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TWA 30B 


TWA 32 




Fig. 1 . Spectral energy distributions for TWA 30B, 32, 33, and 34. Data points are 2MASS (green), ALLWISE (blue), Herschel (purple; see Liu et 
[ak 2015[ l, and ALMA (orange). Two blackbody dust disks are fit in each case. The fractional luminosities, r = L IR /L boh are indicated. The poor fit 
to TWA 30B suggests a more complex model is required. 


2. Observations 


We carried out an ALMA survey of 15 members or proposed 


members of the TWA as drawn from 

Looper et al.) (2007, 

2010# 

l; Looper (J2011); 

Shkolnik et al. 

(j2011); Rodriguez et 

al. (201 

); Schneider et al. 

2012b); Gagne et al. (2014b). These 


stars, listed in Table [1J were chosen as targets as they consti¬ 
tute some of the lowest mass members suggested to date for the 
TWA. None of these had yet been observed with ALMA, though 
most are known to host dusty circumstellar disks as inferred 
via WISE and Herschel infrared excesses (eg, |Schneider et al.| 
|2012a|b[|Liu et al.|2015| l. For those not in published IR surveys, 
we noted the WISE color WT — W4 and marked those with col¬ 
ors redder than 2 as having an IR excess. Some targets, notably 
J1326-5022 and TWA 29, have low TWA membership probabil¬ 
ities (<10%) according to BANYAN II ( |Gagne et al.|2014a l, but 


were nevertheless included in our survey. BANYAN II returns 
a 98% likelihood of membership for TWA 31; however, other 


studies do not consider it a likely member (e.g., Ducourant et al. 

12014) . 


Our ALMA Cycle 2 program (2013.1.00457.S) consisted of 
observations of continuum dust emission at 230 GHz and ob¬ 
servations of the i 2 CO(2-1) and l3 CO(2-l) emission lines with 
a resolution of 488 kHz, corresponding to a velocity resolution 
of 0.6 km/s. We reached a sensitivity of 0.05 mJy/beam in the 
continuum and 5 mJy/beam per 0.6 km/s channel in 12 CO and 
13 CO. Calibration and cleaning was performed by the ALMA 
staff with CASA version 4.2.2. Briggs weighting was used with 
robust=0.5. The final restored beam was, on average, 1.5 x 0.8" 
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and corresponds to a scale of 40-80 AU at the mean distance of 
the stars in the TWA. 


3. Results 


Of the 15 targets observed, four systems were detected in con¬ 
tinuum emission. All four sources were unresolved. These first¬ 
time detections are listed in Table [T] along with the measured 
flux from the primary beam corrected images. Fluxes were mea¬ 
sured by fitting a Gaussian to the cleaned image. We also list 3-<x 
limits, which are three times the RMS error estimated in each 
case. The measured fluxes are consistent with emission from 
T~40 K dust grains as determined by modeling the spectral en¬ 
ergy distribution (SED). Figure[T]shows SEDs of these 4 systems 
along with blackbody fits to their IR/submm excesses (see, e.g., 
Schneider et al.|2012a|. In the case of TWA 30B, which displays 


strong variability and may have complex disk structure as a re¬ 
sult of the edge-on disk geometry ( Looper et al.[2010bj Principe 
et al., submitted), a more sophisticated model is clearly required 
to accurately describe its SED. 


These SED models are overly simplistic; unless the emis¬ 
sion is coming from narrow rings, the real disks will have a 
range of temperatures. Nevertheless, the models in Fig. |T| are 
useful to demonstrate the presence of cold dust in the system. 
The fractional luminosity, r = L/« /L/,,,/, ranges from 2 to 4% 
for TWA 32, 33, and 34. If we assume the continuum emission 
is coming from large (mm-sized) grains which radiate as black- 
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Name 

RA 

Dec. 

Sp. 

Type 

W1-W3 

(mag) 

W1-W4 

(mag) 

IR 

Excess 

Distance 

(pc) 

Flux 

(mJy) 

M cl list 

(10- 2 m £ ) 

Ref. 

TWA 3 0B 

11:32:18 

-30:18:31 

M4 

4.96 

7.23 

Y 

45.8 + 4.8 

0.83 ± 0.07 

3.7 

1 

TWA 30A 

11:32:18 

-30:19:51 

M5 

1.72 

3.64 

Y 

45.8 + 4.8 

<0.20 

<0.9 

2 

TWA 31 

12:07:10 

-32:30:53 

M5 

1.80 

3.92 

Y 

55.0 + 6.4 

<0.20 

<1.3 

3 

TWA 33 

11:39:33 

-30:40:00 

M5 

1.63 

3.25 

Y 

46.6 ± 5.2 

0.33 ± 0.03 

1.5 

4 

TWA 34 

10:28:45 

-28:30:37 

M5 

1.35 

2.75 

Y 

47.0 + 5.6 

0.54 ± 0.06 

2.5 

4 

TWA 32 

12:26:51 

-33:16:12 

M6 

1.77 

3.77 

Y 

59.4 ± 6.2 

2.10 + 0.05 

15.8 

3,5 

J1045-2819 

10:45:52 

-28:19:30 

M6 

0.32 

<2.73 

? 

52.6 ± 6.4 

<0.15 

<0.9 

6 

Jill 1-2655 

11:11:28 

-26:55:02 

M6 

0.47 

0.93 

N 

43.8 + 5.0 

<0.16 

<0.7 

6 

J1203-3821 

12:03:59 

-38:21:40 

M8 

0.80 

<3.09 

7 

65.8 ± 8.6 

<0.15 

<1.4 

6 

J1252-4948 

12:52:09 

-49:48:28 

M8/9 

0.86 

3.00 

Y 

91.8 + 9.4 

<0.13 

<2.3 

6 

J1106-3715 

11:06:44 

-37:15:11 

M9 

<0.88 

<4.02 

? 

53.4 + 8.0 

<0.12 

<0.7 

6 

J1326-5022 

13:26:53 

-50:22:27 

M9 

1.44 

<3.32 

7 

71.0 + 9.0 

<0.13 

<1.4 

6 

J1247-3816 

12:47:44 

-38:16:46 

M9 

2.16 

4.27 

Y 

63.8 + 6.4 

<0.15 

<1.3 

7 

TWA 29 

12:45:14 

-44:29:07 

M9/9.5 

<0.43 

<3.90 

? 

33.3 + 7.2 

<0.14 

<0.3 

8 

J1207-3900 

12:07:48 

-39:00:04 

L0/1 

0.43 

<4.44 

? 

60.2 ± 5.2 

<0.14 

<1.1 

7 


Table 1. Targets for our ALMA observations. The third column indicates if the source had prior indications of IR excess from WISE (WI — W4 > 2); 
"?" denotes upper limits at WISE bands IV3, W 4, or both. Distances listed here have been calculated via the BANYAN II kinematic tool ( |Gagne| 
|et al.|[2014a| l. Also listed are the continuum detections and 3-cr upper limits. Dust masses are estimated as described in Section [3] and using 
Tdust ~ 40 K. All sources are unresolved with ALMA. 

References: (1) Looper et al.|2010b] (2) 


|Looper|201l[ (7) |Gagne et al.|2014b| (8) 


Looper et al. 


Looper et al. 


2010a (3) Shk olnik et al.|20n| (4) |Schneider et al.|2012b| (5) |Rodriguez et al.|2011| (6) 


2007T 


bodies at ~40 K, then the dust grains need to be located a few 
AU from these low-mass stars. 

In addition to the four Table Q] stars that are coincident with 
submm continuum sources, we detected continuum sources that 
are well offset (typically by ~10") from three systems (TWA 
31, J1247-3816, and J1207-3900). Given the typical positional 
accuracy of ALMA observations is better than 0.1", we conclude 
these are background sources, likely of extragalactic origin. 

Assuming optically thin dust, we can estimate the dust mass 
from: 


M (lust — 


F V D Z 


k v B v (T d us t) 


where we adopt k v = 1.15 cm 2 /g (see Rodriguez et al. 2010 and 
references therein) and Td ust = 40 K. We tabulate the resulting 
estimates of Mj, as well as 3-cr upper limits, in Table |T] The 
inferred dust masses range from ~1 to 16% of an Earth mass. 

An alternative approach to determine dust masses is to use 
the temperature-luminosity relationship for protoplanetary disks 


derived in Andrews et al. (2013), namely T r j » 25(L*/L G )^ 4 . 
While consistent with results for disks orbiting young, earlier 
type stars, the |Andrews et al.] ( j201 3j i relationship does not nec¬ 
essarily hold for very low-mass objects such as those in Table [I] 
(van der Plas et al., in prep). Young M5 stars have log L*/L Q ~ 
-2, which suggests a dust temperature of ~8 K. This is 5 times 
lower than the temperatures assumed in the SED model and, un¬ 
der the same assumptions, implies dust masses that are 5 times 
higher. 

In addition to obtaining flux upper limits for the individual 
continuum non-detections, we created a stacked image by aver¬ 
aging non-detections, to assess whether a significant fraction of 
them might display emission at levels just below detectability. 
No detection is apparent in this average image; we obtain a 3-cr 
upper limit of 0.05 mJy/beam. 

Among our sample only TWA 34 displays detectable 12 CO 
emission, which we discuss in Section 0 For those systems 
in our sample with no CO detections, we can infer a 3-cr upper 



Fig. 2. Velocity map for the 12 CO(2-l) emission in TWA 34. 


limit of -0002 M E in molecular gas following the prescription 
in Section 3.1 and assuming optically thin 12 CO, CO:EL of 10~ 4 , 
and a distance of ~60 pc. 


3.1. TWA 34: CO Detection 

Among the Table |T] systems, only TWA 34 shows evidence of 
12 CO emission. Although the l2 CO emission is unresolved for 
each individual velocity channel, we find the centroid changes 
with velocity, allowing us to generate the first moment map 
in Figure [2] Figure [2] suggests that TWA 34 is orbited by a 
molecule-rich disk viewed at intermediate to high inclination, 
with North-South rotation. 
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We show in Figure [3] the integrated line profile of the CO 
emission, Hanning smoothed with a kernel size of 3 channels. 
This double-peaked CO line profile is indicative of Keplerian ro¬ 
tation. Hence, to characterize this emission, we fit a parametrized 
Keplerian model as described in jKastner et al. ( [2008} >. That is, we 
fit a parametric line profile function described by: 


F = 


F 0 ((v-v 0 )/v d ) 3 *- 5 
Fo ((v — vo)/v d ) Pd 


|v - Vo I > V d 
\v - v 0 | < v d 


where F o is the peak line intensity, vo is the rest frequency in 
the star/disk system frame, v d is the projected rotational veloc¬ 
ity near the outer edge of the disk, and p d and q are quasi¬ 
physical disk parameters (see Kastner et al.||2008 for details). 
We fix <7=0.5 and p d - 0.1 for simplicity. The resulting best-fit 
parametric line profile is displayed in Figure[3] We obtain a peak 
intensity of 0.034+0.002 Jy, a systemic velocity of 2.3+0.1 km/s 
in the LSR frame, v d of 2.49 ± 0.09km/s, and an integrated in¬ 
tensity of 0.34 ± 0.03 Jy km/s. The parameter v d can be used to 
estimate the outer radius of the disk as detected in CO emission, 
R d , from 

v 2 d = GM,/R d 


where M, is the mass of the star, 0.08 M a ( jBaraffe et al.|2015) . 
We thereby estimate the CO disk orbiting TWA 34 is ~11 AU 
in radius. This suggests that the disk around TWA 34 is more 
compact than those seen around younger brown dwarfs and low- 
mass stars in Taurus ( |Ricci et al.||2014[ ). We note that the CO 
emission in Figure[2]is marginally resolved and suggests a larger 
CO radius of '-20—40 AU. This discrepancy in the CO outer ra¬ 
dius estimates can also be seen in other disks when comparing 
single-dish line measurements to resolved interferometric imag¬ 
ing (e.g., Kastner et al.|2008| Rodriguez et al.|2010| Sacco et al 


|2014||Huelamo et al.|2015j >7Higher resolution CO imaging will 

allow an independent measurement of the disk size, which can 
in turn be used to accurately estimate the mass of TWA 34. 

We estimate a gas mass of ~0.2 M E following the pre¬ 
scriptions of Zuckerman et al. (2008), Kastner et al. (2008), 
Rodriguez et al.l (|2010)), and references therein. We assume 
optically thin _T3 C071T ,2 C: ,3 C ratio of 89, a CO:H 2 ratio of 
10 4 , temperature of ~40 K, and a 3-cr 13 CO upper limit of 
~0.03 Jy km/s. The 13 CO upper limit is determined assuming 
a linewidth identical to that of 12 CO. This estimate for the gas 
mass is ~1 times higher than that for the dust mass. This ra¬ 
tio is comparable to other evolved molecular gas disks such as 
TW Hya and V4046 Sgr, as well as some disks in the much 
younger (~l-3 Myr) Taurus star-forming region (e.g., |Sacco et 
|al.|2014j Williams & Best [2014] ). As these and many other previ¬ 
ous studies suggest, such low inferred gas-to-dust ratios could be 
indicative of gas removal through accretion or photoevaporation, 
or could merely reflect overestimates of the CO: Hi ratio; due, for 
example, to freezing out of CO on cold dust grains. These var¬ 
ious processes result in large uncertainties when estimating the 
gas mass of disks based on CO measurements. 


3.2. TWA 34: System Velocity 

From our best-fit Keplerian profile (see prior section), we have 
obtained an estimate for the systemic velocity of TWA 34. In 
the barycentric frame of reference, this velocity corresponds to 
13.3 ± 0.1 km/s. This is the first accurate radial velocity mea¬ 
surements available for TWA 34 (see also [Murphy et al.]2015| l. 
At a distance of 47.0 + 5.6 pc and using the proper motions listed 
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Fig. 3. 12 CO(2-l) emission line profile of TWA 34, the only target in 
our sample with detected CO gas. The red, thin line represents the best- 
fit Keplerian profile. 


in PPMXL ( [Roeser et al.|201Q] >, we can estimate UVW veloci¬ 
ties of -11.0 ± 1.6, -16.2 ± 0.6, -3.9 ± 1.4 km/s. These agree 
very well with the average velocity of the TWA (-9.87, -18.06, 
-4.52 km/s; Malo et al.|2013| ). This new velocity measurement 
thereby further supports the conclusion that TWA 34 is a mem¬ 
ber of the TWA. 


4. Discussion 


Among the 15 low-mass TWA members and candidate mem¬ 
bers listed in Table [T| only 4 yielded ALMA continuum detec¬ 
tion at 1.3 mm, despite the fact that many show some evidence 
of warm circumstellar dust. In the absence of cold dust grains, 
the warm dust grains detected by WISE would have 1.3 mm 
emission <0.1 mJy, which is below the sensitivity of our ALMA 
observations. Because we are only sensitive to cold mm-sized 
grains, our observations appear to demonstrate that the presence 
of warm circumstellar dust does not necessarily imply cold dust 
is also present in the system. This is in agreement with prior 
studies of M stars (see, e.g., Lestrade et al.|2009] and references 
therein) that have found the incidence of cold disks to be much 
lower around such stars than around higher mass stars. Although 
the ALMA non-detections appear to indicate no cold grains ex¬ 
ist, an alternative explanation is that any surviving grains in the 
outer disk have already grown to cm size or larger and become 
invisible at 1.3 mm wavelengths (Ricci et al.|2010a|b[|Mohanty 
let al.|2()13| l. 

The continuum non-detections suggest dust masses of about 
10 ~ 2 Me or less. This is similar to what has been observed for 
other ~ 10 Myr-old debris disks around M stars (see |Wyatt|2008| 
and references therein). We were more sensitive to molecular 
gas masses and achieved a limit of a few times 10 3 Me for gas 
in H 2 , assuming CO/H 2 of 10 4 . As the case for TWA 34 and 
other disks shows, the gas-to-dust ratio is unlikely to be ~ 100 
as in the ISM (see Williams & Best |2014] l. It appears likely that 
by the age of the TWA, the gas in a typical M star’s disk has in 
general been efficiently removed, even in cases where a signifi¬ 
cant mass of primordial dust has survived. However, studies have 
identified signatures of on-going gas accretion in some of these 


systems, for example around TWA 30A, 30B, and 31 (Looper 


|et al.||2010a|b| [Shkolnik et al.|[20TT) . Hence, at least in certain 
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cases, it is likely that disk CO gas has frozen out onto dust grains, 
suppressing the gas-phase CO abundance. 

5. Conclusions 

We have carried out an ALMA survey of 15 low-mass TWA 
members and candidates to search for molecular gas in the form 
of 12 CO and 13 CO as well as provide constraints on continuum 
dust emission. Among systems targeted, four (TWA 30B, 32, 33, 
and 34) have detected dust emission consistent with the existence 
of cold dust grains in the disk. Circumstellar dust grain temper¬ 
atures of ~40 K are consistent with the mid-infrared to submm 
SEDs for these systems. All continuum sources are unresolved. 
While most of our sample shows indications of warm dust based 
on WISE measurements, the ALMA non-detections suggest any 
cold grains present in the outer disk may have already grown to 
cm size or larger. 

Only one system, TWA 34, shows signatures of molecular 
gas in its disk in the form of 12 CO (2-1) emission. The 12 CO 
emission has velocity structure indicative of Keplerian rotation. 
The systemic velocity for the system, as determined from the CO 
detection, is consistent with membership in the TWA. Among 
the sample of known ~7-10 Myr-old star/disk systems, TWA 34, 
at just ~50 pc from Earth, is the lowest mass star thus far identi¬ 
fied as harboring cold molecular gas in an orbiting disk. 
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